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Abstract 

The Non-thermal phase transition in high energy colli- 
sions is studied in some detail in the framework of random cas- 
cade model. The relation between the characteristic parameter 
Xq of phase transition and the rank q of moment is obtained 
using Monte Carlo simulation, and the existence of two phases 
in self-similarly cascading multiparticle systems is shown. The 
relation between the critical point Qc of phase transition on the 
fluctuation parameter a is obtained and compared with the ex- 
perimental results from NA22. The same study is carried out 
also by analytical calculation under central limit approxima- 
tion. The range of validity of the central limit approximation 
is discussed. 
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Recently, the prediction'^1 that there exist the property of self-afiine fractal in the 
anisotropic phase space of multiparticle final states in high energe hadron-hadron colli- 
sions has been confirmed by experiments^^'^'. This breakthrough in the nonlinear study 
of high energy physics places the further study of nonlinear property of multiparticle 
final states on the agenda. 

In this respect, the non-thermal phase transition l^'^l is a prolem worthy while fur- 
ther study. In the presently available experiments'^' , due to the restriction of energy, 
the average multiplicity is very low, and the rank of the factorial moments could not 
be high. So, no clear evidence of non-thermal phase transition has been seen. The new 
Large Hadron Collider (LHC), which is being built and will be put into operation in 
the beginning of next century, will dramatically raise the collision energy and multi- 
plicity, providing perfect condition for the study of non-thermal phase transition. For 
a theoretical preparation it is necessary to carry on detailed discussion on this phase 
transition and to clarify its property. 

The aim of this short paper is to make a model study of the non-thermal phase 
transition, especially to make clear of the relation between the critical point of non- 
thermal phase transition and the strength of dynamical fiuctuations. 

The random cascading a model is widely used in the study of nonlinear property 
of multiparticle final states in high energy collisions. Using this model, it is easy to 
get a system pocessing the property of intermittency and fractal. We will show that 
non-thermal phase transition does exist in this system and the relation between the 
critical point of phase transition and the parameter a of fiuctuation strength in the 
model can thus be obtained and compared with the experimental data. 

Firstly, let us briefly remind the random cascading a model'^l with probability 
conservation. 

Consider a region A of one-dimensional phase space. Devide it into A cells. The 
probability of particles faUing into the ith cell is 

Pi^PQUJi, (1) 

where Pq = 1 is the probability in the phase space region A, uji is the probability of 
the elementary partition. Next, we divide each sub-bin into A even smaller sub-bins. 
The probability in the ijth bin (i = 1, 2, . . . , A; j = 1, 2, . . . , A) is 

Pij^Pii^j, (2) 
After u steps, the probability in a sub-bin is 

V 

Pili2-i. = Il^tk- (3) 

fe=i 



The total number of intervals is M = . 

In order to guarantee the conservation of probability in each step of cascading, we 
choose the elementary probability uji for A = 2 asl^^: 

1 + ar 1 — ar 

^1 = — > ^2 = — . (4) 

where is a random number distributed uniformly in the interval [—1, 1], a is a model 
parameter describing the strength of nonlinear dynamical fluctuations (0 < a < 1). 
The definitions of the probability moments and factorial moments are'^' : 

c, = f^^;;=^^-^ = m^-\j:pI). (5) 

[M\^m=lPm}) m=l 

It can easily be proved that under the assumption of Poisson or Bernoulli type of 
statistical fluctuations the normalized factorial moments Fg are equal to the normalized 
probability moments Cg. 

The character of dynamical fluctuations can be expressed as the anomalous scaling 
of probability (or factorial) moments: 

Cg{M) oc M^^ 

or equivalently 

\nCg(M) ^ A + ipg\nM (M ^ oo) . (7) 

where (fg is called intermittency index. 

In order to see the anomalous scaling of probability moments more clearly, we 
choose the fluctuation-strength parameter a — 0.5, the elementary partition number 
A = 2, the division step u — 12, the ranks of moment g = 5, 10, 15, 20, 25, 30, and make 
use of Eq.(5) to simulate the relation lnCg~lnM. The results are shown in Fig.l. 
The intermittency parameters (pg are obtained through linear fit. We can see from the 
figure that the higher the rank q is, the larger the slope (pg is. 

A parameter A^ has been introduced '^'^1 in the multifractal analysis to characterise 
the non-thermal phase transition in the multiparticle systems. It is related to the 
intermittency index (pg by the relation 

Xg = {(Pg+l)/q. (8) 



We will try to evaluate this parameter both through analytic calculation and by using 
Monte Carlo simulation. 



In the random cascading a model, the probabihty moment is: 



\^) 

It can be rewritten as: 

= X'^''(exp(^-j2qe)jy (10) 

where Ei — —Incu (i). The parameter ( — Yli=i Q^i in the above equation is the sum of 
u random numbers. Under the central limit approximation ( approades to Gaussian 
distribution: 

C, (M) = A^''(e-^^) = exp ^qu In A + ^ - ■ (H) 
Using Ci(M) = 1, we get 

Cg (M) = e^'^'''(«-i)/2_ ^^2) 
The intermittency indices can be deduced as: 

<P, - ii^il^. (13) 
21n2 ^ ^ 

We have also the relation 

= {In^u) - {Inuf ^\a'' ^\a^ ^ . 

Under linear approximation'^^ it becomes cr^ = Substituting into Eq.(13) we get 

= iii^l^. (14) 
61n2 ^ ^ 

from eq.(8): 

(^, + 1 (g-lK 1 
g In 2 g 

The resulting Ag~g are plotted in Fig. 2(a) for a = 0.2,0.3,0.4,0.5 respectively. 

Fig. 2 (a) is the result under central limit approximation. The exact relation can not 
be calculated analytically. Therfore, we use Monte Carlo simulation. The resulting 
Ag~g for a = 0.2,0.3,0.4,0.5 respectively, are shown in Fig. 2(6). 

It can be seen from the figures that the Xq^q curves from both the Monte Carlo 
simulation and the analytical calculation under central limit approximation have the 
same trend, i.e. with the increasing of q, \g arrive at a minimum at the point qc, 
which means that there really exists non-thermal phase transition in the self-similar 
cascading model and two different phases do indeed coexist, qc is the critical point of 
phase transition. 



In Fig. 2 we alse draw the experimental data from NA22 (open circles). It stops 
at the rank g = 5 and is unclear whether there is a minimum at some higher rank as 
required by non-thermal phase transition. The open triangles in the figure is the result 
from the same experiment selecting only the particles with low transverse momenta 
{pt < 0.15 GeV/c). In this case, with the increasing of q (from 4 to 5), increases. 
It seems to show that there is phase transition and the critical point Qc < 5. As is 
well known, choosing only the particles with low transverse momenta, the strength of 
intermittency increases'^*^' . Therefore, this experimental phenomenon shows that the 
system with lower transverse momenta, which has larger intermittency strength, has 
lower critical point of non-thermal phase transition. This is qualitatively the same as 
the result of our model, where the phase transition point shifts left with the increasing 
of fluctuation strength (when a increases Qc decreases). 

In order to see more clearly the relation between the phase transition point Qc and 
the fluctuation parameter a of the model, we draw the figure of Qc'^a, as shown in 
Fig. 3. We can see from the figure that the larger a is, the earlier Qc appears. 

Comparing the exact values of gc~CK from Monte Carlo simulation and the analyt- 
ical results under central limit approximation, it can be seen that both have the same 
trend of continuously descending. However, the values of in central limit approxi- 
mation are generally smaller than the exact values. This shows that the central limit 
approximation can refiect qualitatively the property of non-thermal phase transition 
but there is noticeable quantitative deviation. 

References 

[1] Wu Yuanfang and Liu Lianshou. Phys. Rev. Lett., 1993, 70: 319; Science in China, 
1995, A38: 4357 

[2] EHS/NA22 CoU., Agababyan N M et al. Phys. Lett., 1996 B382: 305; ibid., 1998 
B431 451 

[3] Wang Shaoahun, Wang Zhaomin and Wu Chong. Phys. Lett., 1997 B410: 323 
[4] Peschanski R. Nucl. Phys., 1989 B327: 144 
[5] Bialas A, Zalewski K. Phys. lett., 1990 B238: 413 
[6] De Wolf E A, Dremin I M and Kittel W. Phys. Rep., 1996 270: 1 
[7] Wu Yuanfang, Zhang Kunshi and Liu Lianshou. Chinese Science Bulletin, 1991 36: 
1077 

[8] Bialas A and Peschanski R. Nucl. Phys., 1986 B273 703; ibid., 1988 B308 857 
[9] Liu Lianshou, Fu Jinghua and Wu Yuanfang. Phys. Lett., 1998 B444: 563 
[10] Wu Yuanfangand Liu Lianshou. Phys. Lett., 1991 B269: 28 



o 

c 



80 
70 
60 
50 
40 
30 
20 
10 




¥ q = 30 
A q = 25 
• q = 20 
□ q=15 
A q=10 
q = 5 



T 
A 
• 

□ 
▲ 



T 

A 







T 


A 


• 


A 


• 




• 


□ 


□ 


□ 
▲ 


▲ 


▲ 






1 




1 


2 


3 





A 

□ 



A 



□ 



▲ 





T 

A 



□ 



A 



□ 



A 



□ 







A 



□ 



A 



□ 







4 



5 6 



9 

nM 



Fig.l Log-log plot of various rank probability moments versus partition number in a 
model 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.4 

0.35 
0.3 

0.25 
0.2 

0.15 
0.1 



a = 0.5 
a = QA 
a = 0,3 
a = 0.2 



a 




b 



a = 0.5 



a = 0.4 



a = 0.3 



T T T T ' ' 



c( = 0.2 



5 10 15 20 25 

q 



Fig. 2 Relation between the parameter and the rank q of moments. The vertical 
lines indicate the position of minima, (a) Analytical results under central limit ap- 
proximation. The open circles are the experimental results from NA22. Open triangles 
are the results from the same experiments taking only low momentum particles with 
{pt < 0.15 GeV/c) Data taken from Ref.[6]. (6) Results of Monte Carlo simulation. 
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3 The relation between phase trnasition point Qc and fluctuation strength 



